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A formal solution is obtained for the magnetic field produced by a vertical magneto- 
static dipole which moves in free space at a constant velocity along the normal to a con- 
ducting half space. When the velocity of the dipole is much less than that of light, the 
field in free space is governed by Laplace's equation, and that in the conducting half space 
by the diffusion equation. An expression in closed form is derived for the voltage induced 
in a small loop moving with the dipole. For certain values of the parameters, simpHfied 
expressions are obtained for the induction. Numerical results are presented for certain cases. 



1. Introduction 

The problem of a. stationary vertical oscillatinsi- 
magnetic dipole above a conducting ground has been 
discussed by Wait [1951, 1953]. The analagous 
problem of a vertical magnetostatic dipole approach- 
ing a conducting liall' s])ace is considered in the 
present paper. If tbe dipole moves with a constant 
velocity which is much smaller than that of light, 
then the magnetic field in free space will satisfy 
Laplace's equation while the fiehl in the con- 
ducting region will satisfy the diffusion equation. 
Formal solutions are obtained for the magnetic 
field in the two regions. 

It is of interest to find an expression for the vol- 
tage induced hi a small loop moving with the dipole. 
Such a result might have application in the design 
of a device for measuring the close proximity of a 
vehicle to a conducting surface or for measuring 
the velocity of a vehicle just before impact with 
such a surface. The solution for the reflected 
vertical component of the magnetic flux density 
in free space is used to find an expression in closed 
form for the induction. This result simplifies 
considerably for certain values of the parameters. 
Numerical results are presented for certain cases. 

2. Formal Solution 

With reference to the coordinate system shown 
in figure 1, the conducting medium of conductivity 
0-, permeability /x, and permittivity e occupies the 
half space defined by all negative values of z. The 
upper half space has the ciiaracteristics of free 
space, namely, zero conductivity, permeability /zq, 
and permittivity €o. A magnetostatic dipole of 
moment m moves toward the conducting half 
space with a constant velocity v along the z axis. 
The dipole is situated instantaneously on the z 
axis at z=d and is oriented in the positive z di- 
rection. The field point P, located at p, <^, z, in 



cylindrical coordinates, is a distance r from the 
dipole. The vector magnetic flux density B can 
be represented everywhere in terms of a magnetic 
vector potential A. 
Thus, 



P=VX^ 



(1) 



wliere /I, if i <Cc (c = free space velocity of light), 
satisfies 



V'vl=0 



for 2>0, and 



for 2<0. 



VM= 



= M0- 



dt 



(2) 



(3) 
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Figure 1. Magnetostatic dipole moving toward a conducting 
half space. 
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In the region 3>0, the total vector potential is 
assumed to consistjof a primary potential Ao plus a 
reflected potential Ai. The primary vector potential 
of the magnetostatic dipole, shown in figure 1, has 
only a component which is given by 



-m 



Mo 



47r dp 



(4) 
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The quantity - may be expressed [Smythe, 1950] as 



where r= 
then reads 

Ad,~-=m 



^ Jo 
[p^-\-(d—zyY^^. The primary potential 



Mo [ 

47rJo 



e-'^'-'^J,{kp)kdk, (3>0). (6) 



Similarly, the reflected potential in the upper lialf 
space, ^0p and the potential in the conducting half 
space, ^02? may be written as 



A^^=m 



Mo 
47r 



r i,{k)e-'^'+''J,{kp)kdk, 
Jo 



(2>0) 



and 
A 
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= m -r- 
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e-fa^e-'"'Ji(kp)k(lk, 



(2<0), 

(8) 
respectively. 

The functions a{k),fi (k), and/2(/:) are determined 
from (2) and tlie boundary conditions at 3=0 for the 
normal magnetic flux density and the tangential 
magnetic intensity. 

Upon substituting the expression for A^^ into (3), 
the relationship 

a^=^k^+k^av (9) 



is obtained where v, the velocity of the dipole, 
Application of the continuity condi- 



, dd 
equals-^^ 



tions for the magnetic field at the boundary 2=0 
yields 



Acj)Q-{-A^^ — ^ 



02' 



(normal B continuous) 



and 



(10) 



Mo Mo 



A 



4>2 



(tangential H continuous) . 



(11) 



The functions ji{k) and f2{k) are evaluated by 
substituting the potentials given by (6), (7), and (8) 
into (10) and (11). This results in 



my- 



pk- 



-Mo« 



pk+p^a 



(12) 



and 



Uk) = 



2pk 



pk + poa 



(13) 



The components of the magnetic flux density in 
each region are then found from 



D 1 ^ / ^ \ 

p Op 



and 



B, 



" dz' 



(14) 
(15) 



3. Voltage Induced in Loop Moving 
With Dipole 

3.1. Solution in Closed Form 

An interesting application of the solution for the 
reflected magnetic flux density in free space is the 
determination of the voltage induced in a small, 
single-tm^n loop moving with the dipole along the 
z axis. If the axis of the loop is vertical and the 
permeability of both regions is taken as pq, the 
voltage per unit loop area, V, is given by 



V- 



dB,^ 



mvpo 



dt 



TT Jo 



-(M 



1+ 



(?^' 



Vdk, 



(16) 



where p^=poav. The integral in (16) is evaluated 
in appendix 1 and can be written as 



F- 



mvpoF 3 3 



-2A^(-4^+,' 



where h=2phL 

3.2. Solution for h<d 

For h<^\, the expression for the voltage given by 
(17) is not useful for numerical computation due to 
the fact that the sum of the first two terms is op- 
posite in sign and very nearly equal to the sum of 
the third and fourth terms. A more convenient 
expression may be derived by substituting h=2jM 

k . . 

and ^ = -^ in the integral \\\ (16). The equation 

then reads 



V=mv 



[j; 



p^e ^^^dx 



-X 



+ 2pH^{x-[x'+xY')e-'''dx 



]• 



(18) 
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Evaluation of tlu^ first iritc^o-ral in (18) is simple and 

yields ~- jv^- The second iutepral is evaluated 

loTT a 

in appendix 2 hy neglecting terms of tlie order of Iv" 
and higher. The expression for the voltage then 
simplifies lo 



V= 



mviJiop^_maiJilv^ 



327rd' S2t(P 
3.3. Spacial Case for a= oo 



(19) 



For the case wliere the conductivity of the lower 
half space becomes infinite, tfie solution for the vol- 
tage may be obtained by either tlie method of images 
or b}^ using (16). The resulting expression is 



XT' *^ '^ / \ 

^^i^ ^^MoT4^ (o-=oo). 
1 DTT (r 



4. Discussion of Numerical Results 



(20) 



When /i<l, or equivalently wlien aV(l<^ % X 10^ 
the expression for the voltage given by (19) may be 
used. Numerical results using this expression aie 
presented in figures 2 and )^ for the case where 7n=l 
amp/m^. This value of the moment corresponds 
to that of 1 cni'^ of Alinco 5 in which there is an 
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Figure 2. Induced voltage as a function of height above a 
conducting half space. 



average flux density of 1.2 webers/m^. In figure 2, 
the voltage induced in a loop of unit area is plotted 
as a function of the lieight of the inngnet above 
surfaces of different conductivity. The velocity of 
the loop and magnet is taken as 300 m/sec. Figure 
3 shows the voltage per unit area as a function of 
velocity for various half space conductivities. The 
values of induction in this figure correspond to a 
one meter height of the loop and magm^t above 
the conducting suiface. 

The numerical results plotted hi figures 2 and 3 
show the induction to be quite small. Usable volt- 
ages could be obtained just above earth or sea 
water by employing large magnets and small multi- 
turn loops moving witli velocities near Mach 1. 

5. Appendix 1 

k 
Usuig tlie substitution x=-7^} the hitegral in (16) 

is written as 



V=7nv 



Mo 



' [ [2/e-2-^V-2//(j2+:^)i/2^-2xp2rf^3 

+p'e-'''^^'j'\(Ix. (Al) 
The first and third terms are readilv inteu'rated. 
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Thus, 



mvpLo 



The second integral can be written as 

TT Jo 

_nm^d^r^^,_^^yj,^_n.^^^ (A3) 

where h = 2p^d. If the substitution x=^ cosh /3— ^ 

is made m the integral on the right side of (A3), 
the integral reads 



/^ 00 1 /^ °° 

{x'^ + xy^^e-'''dx=- sinh^iSg 



Hcosh/S-l) 



f/^ 



(A4) 

Since Kn{u), the modified Bessel function of the 
second kind, is defined [Erdelyi, ed., 1953] by 

J» 00 
-"^^^^^coshni86?/3, (u>0), 
0^ 

(A4) can be written as 

j^^ {x'+xy^'e-'^dx=le^ ^K2 (|)-^o Q)] 

The third derivative with respect to h of the integral 
in (A5) is obtained by using the following recurrence 
relations: 

K',{u) = -K,{u) 

where the primed quantities represent derivatives 
with respect to u. The second integral of (Al) is 
finally expressed as 



— mvnQp 



'^^ rx^(x'+xy^h-''dx 

TT Jo 



^niVMt ^1 1 



+(i;-i+2i=-i)^.(l)]- <-) 



6. Appendix 2 

By letting t=hx, the second integral in (18) can be 
written as 



^ f ^x'[x-{x'+xy^']e-'''dx 

TT J 

( ^t'[t-(t'+thy^'e-'dt. (A7) 



r 
mvjjLop' 



The integral on the riglit of (A7) is written in two 
parts, namely, 



I,= f \'[t-(t'+thy'']e-'dt 



and 



I^=r t\t-{t'+thy"]e-'dt. 

The |/i| is at most of the order h^ since 

|/i|=-| {\M-{t'+thy"]e-'dt\<{\\t-2"'h)dt 
Jo Jo 

h' 2'^'h' 



(A8) 



For^<l, 



(t'+thy^'=t{i+j) 

where < 9(0 < 1 . It follows that 

- r T-^^ ^+^ ^%^"| e-'dt-B, (AlO) 

where iZ=T^ ^[l+^(^)] ^^^6 ^dt represents a re- 
io J h 

mainder term. The maximum value of the 

remainder is found by letting d(t) =Q. Thus, 

Therefore, 

(A12) 



Ix+l2=-^hVih?+(i{h^). 
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The second integral in (18) is then written as 

mviJLo 



27r 



bi+ife]- ^^'''^ 



The author thanks Dr. A. H. Van Tiiyl of the U.S. 
Naval Ordnance Laboratory, White Oak, Maryland 
for his help in evaluating the integrals considered in 
the appendices. 
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